ON FUNCTIONS WITH A CONJUGATE 
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Abstract. Harmonic functions of two variables are exactly those 
that admit a conjugate, namely a function whose gradient has the 
same length and is everywhere orthogonal to the gradient of the 
original function. We show that there are also partial differential 
equations controlling the functions of three variables that admit a 
conjugate. 



1. Introduction 

A pair of smooth real-valued functions / and g on a Riemannian 
manifold M are said to be conjugate if and only if 

(1) ||V/|| = \\Vg\\ and (V/, Vg) = 0. 

In this article, we shall address the following question. When does a 
given smooth function / : M — > R admit a conjugate function? When 
M is 2-dimensional the pair of functions (/, g) : M — > R 2 is mutually 
conjugate if and only if the mapping (/, g) is conformal away from 
isolated points where its differential vanishes. It is well-known that, 
in this case, / must be harmonic and, conversely, a harmonic function 
locally always admits a conjugate, unique up to an additive constant. 
When M is of higher dimension, then the pair (/, g) : M — > R 2 is 
said to be semiconformal. As discussed in [Bj, semiconformality is one 
of the two conditions that (f,g) be a harmonic morphism. In fact, 
if M = R n and both / and g are polynomial, then it is the only 
condition [!]• In this article, we shall be concerned with / defined on an 
open subset in R 3 . We extend our earlier work [2] in which we derived 
some necessary conditions on / in order that it admit a conjugate. 
Under some mild non-degeneracy assumptions, we now derive necessary 
and sufficient conditions. 
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An example of a pair of conjugate functions in three variables is 

, xi 2 + x 2 + x 2 x 1 2 + x 2 2 + x 3 2 

f = x 2 y~ — 2 — g = x 3 — — - 2 — . 

x 2 + x 3 z x 2 + £3 

The Hopf mapping S 3 — > S 2 viewed in stereographic coordinates 
(1- ||a;|| 2 )x2 + 2xix 3 _ (1 - ||x|| 2 )x 3 - 2x x x 2 

x 2 2 + x 3 2 x 2 2 + x 3 2 

provides another good example. In these two cases, the pair (/, g) 
enjoys an evident symmetry with respect to rotation about the xi-axis. 
This is not usual, as is illustrated by the following example:- 



/ = log y/ Xi 2 + x 2 2 + x 3 2 g 



arccos ■ 



y/xi 2 + x 2 2 + x 3 3 

In all three examples, the pair (/, g) is smooth away from the Xi-axis. 

We shall frequently need to manipulate tensors and for this purpose, 
we use Penrose's abstract index notation |13|. We shall write 



fi = Vi/ fij = ViVjf et cetera, 

where Vj is the flat connection on IR n or, more generally, the metric 
connection on a Riemannian manifold. Also, let us 'raise and lower' 
indices with the metric in the usual fashion and write a repeated index 
to denote the invariant contraction over that index. Thus, /\ = Af 
is the Laplacian and f l gi = (V/, Vp). We shall use round and square 
brackets to denote symmetrising and skewing over the indices they 
enclose. For example, (fi(ij)k = \<fiijk + \4>jik an d ^[i4>j] is the exterior 
derivative of a 1-form 0j. 

2. A NECESSARY CONDITION 

Theorem 1. Let M be an 3- dimensional Riemannian manifold and 
f : M — > R a smooth function. In order to admit a conjugate, f must 
satisfy the differential inequality 

(2) Vffjffh - f l hf jk f Jk + ffiiFtf < 0. 

Proof. A proof of this theorem was given in [2]. In fact, a version was 
proved there valid in any dimension. Here we give a more efficient 
proof only valid in three dimensions. However, this proof will allow us 
to draw additional and useful conclusions. In addition, the method of 
proof will provide a good illustration of the normalisation techniques 
occurring throughout the rest of this article. 

If / is to admit a conjugate, then there must be a closed 1-form Uj 
so that 

(3) f J Uj = and u J ujj = f 3 fj. 
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Indeed, ([I]) implies that we may find an ojj that is exact. We shall 
show that the inequality (|2J) is necessary in order to find a closed un- 
satisfying (J3]). To proceed, let us differentiate the equations with V 1 . 
We obtain 



(4) f^ + ^/, = and J'-'Ij- 

Since we are supposing that u>y = ViUj is symmetric we may transvect 
the second of these with /, and use the first to eliminate u/^/j. This 
gives 

+ r/</i = o. 

We now have the following equations 

(5) fu t = uM = ./•'/, + rVi/i = o 

and we claim it is a matter of algebra to show that the inequality (J2j) 
must hold if there is to be a solution u>j. This follows immediately from 
the following Lemma. □ 

Lemma 2. If is a 3x3 symmetric matrix and is a 3-vector, then 

(6) ffMfiffk ~ ffif jk fjk + ffiiPif] + l2T ijk T« k = 
where 

(7) Tijk = f[iUjfk]eUJ £ 

and u>i is any solution, real or complex, of the equations ((Sj) . 

Proof. The quantity on the left hand side of (|2J) and in the statement of 
this lemma will arise repeatedly in this article. It is useful to introduce 
the notation X for it. If fi = then the conclusion is trivial. Otherwise, 
let us choose coordinates so that f\ = fi = 0. We may also orthogonally 
diagonalise the quadratic form restricted to the plane orthogonal 
to /j. In other words, we may further change coordinates to arrange 
that f\2 = 0. Having made these choices, the quantity X becomes, 
after a short calculation, 

(8) X = 2(/ 3 ) 2 (/ll + /33)(/22 + /33). 

Another short calculation yields 

(9) T l]k T^ k = i(/ 22 - / n ) V V 
whilst the equations (JSJ) become 

(10) u; 3 = Wl 2 + u; 2 2 = /3 2 fiW + hW + / 33 / 3 2 = 
the second two of which may be written as 

fin 1 1 

fll + /33 fl2 + ^33 









\ fs 2 ] 




u 2 2 
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Now there are two cases. If /n = fa, then ((9]) implies that T^T 1 ^ = 0. 
But ( ITT]) implies that /n + ^33 = and then (jSJ) shows that X = and 
()6]) reduces to = 0. On the other hand, if fn 7^ ^22, then we may use 
(fTT]) to solve ( TTOj) . obtaining 



/ 10 \ 2 <■ 2 /22 + /33 j 2 j.2/ll + /33 

(12) wi = / 3 t 7- and w 2 = / 3 -7 5-. 

J22 — J11 J11 — J22 

and compute 

l2T ljk Ti k = 2(/ 22 - / u ) W = -2/ 3 4 (/n + /33X/22 + / 33 ). 
A comparison with (jSj) immediately yields (jBJ), as required. □ 

From now on we shall suppose that fi is non-zero (at a particular 
point and hence nearby as well). In case that / admit a conjugate, 
it is then clear from ([1]) that the pair (f,g) is a submersion (near the 
point in question). The nature of the singularities of a semiconformal 
mapping is not known in general [3J. 

Notice that it follows from the proof of this lemma that the equations 
§5§ always have solutions if we allow u>i to be complex and generically 
(in fact, precisely when 1^0) there are four solutions. Alternatively, 
this is geometrically clear: the first equation restricts matters to a plane 
wherein the second and third equations describe planar quadrics. 

Perhaps our proof of Lemma [2] seems bizarre but, in fact, we have 
used a familiar technique. The Cayley-Hamilton Theorem, for example, 
is often proved, even for real matrices, by employing Jordan canonical 
form over the complex numbers. Not only that, but Lemma [2] can 
be proved without normalisation by means of the Cayley-Hamilton 
Theorem applied to restricted, as a symmetric form, to the plane 
orthogonal to fi (the details of this proof being left to the reader). 
Another proof avoiding normalisation may be obtained by expanding 
the identity = f\iUjfk k fef- In fact, it is a consequence of Weyl's 
Second Fundamental Theorem of Invariant Theory [15] that dimension- 
dependent identities must arise by 'skewing over too many indices'. To 
use normalisation as we have done, however, is a simple enough method 
that we shall employ throughout this article. 

The quantities occurring in the proof of Lemma [2] suggest other com- 
binations of derivatives with geometric significance. The operator 

(13) / • > z f j fifj + ffifj, 

for example is, up to a multiple, the well-known 3-Laplacian [7] [10] and 
in normal coordinates 

(14) h = h = fu = 
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at a point becomes 

(15) Z = f 3 2 (f u + f 22 + 2f 33 ). 
Also, the quantity J = f l /, is f 3 2 . Therefore, from ([8]), 

Y = Z 2 -2JX 

(16) = /3 4 (/ll + /22 + 2/33) 2 - 4(/ 3 ) 4 (/ll + /33)(/22 + /33) 
= /3 4 (/ll — f^lf 

and we recognise that the vanishing of this expression when X = is 
exactly the criterion discovered in the proof of Lemma [2] for there to 
be infinitely many solutions uoi to the system (151). In summary, if we 
allow complex solutions of fl5]) then 

1^0 3 4 distinct solutions 

(17) X = and Y ^ 3 2 distinct solutions 
X = and Y = -<=>- 3 oo-many solutions. 

If we restrict attention to the case when (jSJ) has real solutions, then 
Lemma [2] implies that X < whence 

1^0 -<=>- X < and 3 4 distinct solutions 

(18) X = and Y ^ F > and 3 2 distinct solutions 

Y = <^=^ X = and 3 oo-many solutions, 

the last two conclusions following from Y = Z 2 — 2JX upon noting 
that both terms on the right hand side are non-negative. 

3. Integrability of the conjugate direction: 
the generic case 

Recall that if / is to admit a conjugate function near any particular 
point, then there must be a solution Uj at that point of the algebraic 
equations (jSJ). These three equations, specifically the third one, were 
derived under the assumption that Uj extend to a closed form near 
the point in question but our approach from now on is to take Uj to 
be defined at a particular point by the equations (jSJ) and ask whether 
it may be extended to a smooth closed form near that point whilst 
maintaining (JS]). This is entirely equivalent to finding a local conjugate 
for /. As a matter of terminology, we shall refer to a solution u>j of 
as a conjugate direction. In case that X < (at the point in 
question and hence nearby as well), we have just seen from (fl8|) that 
there are four distinct solutions of (jSJ) for Uj. It follows that any one of 
these solutions uniquely and smoothly extends as a conjugate direction. 
Therefore, the only remaining question in case X < is whether this 
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extension is closed and we shall refer to this as integrability. We show 
that integrability is equivalent to a further two polynomial equations 
in bj % and the derivatives of /. 

Resolution of these further equations combined with (jSJ) will lead 
to necessary and sufficient differential conditions on the function / in 
order that it admit a conjugate. All of this is under the assumption 
that X < and we shall refer to this as the generic case. The case 
X = will be studied separately. 

Theorem 3. Let Uj be a conjugate direction determined by (jSJ). Then 
provided X < 0, the tensor field Uij is symmetric in its indices if and 
only if 

(19) r k f l f 3 f k + r k f t co J u k + 2Pf J k f l f k -2Pf J k u l u k = 

(20) f^fifjUk + ^UiUjUk + ^tffiUk = 0. 
Proof. Since X ^ 0, the identity of Lemma [2j namely 

(21) /•7,,Y • 127; ; ,7" A - (I. 

where %j k = f^ujjf^uj 1 , shows that the vector field f^coj is independent 
of f % and u;\ Therefore, the tensor field is symmetric in its indices 
if and only if 

(22) u l v j (uij - Uji) = , 

where u l and f- 7 are any vector fields taken from the set {f l ,u\ f^Uj}. 
Looking back at (jlj), which was obtained by differentiating Q, we see 
that 

fuJ^UJij - Up) = f' j f,fj + f ij UiUj. 

This already vanishes by assumption. It is our third equation from (J5]). 
Differentiating this third equation gives 

o = .rvAr'fjh ■ /'^>v) 

= f ijk fifjfk + f jk fiUjU h + + V ik u*fuk ■ 

We notice that the last term fi k Uijf l uJk occurs as the first component 
of the symmetry condition /' p k uJk(^ij—^ji) — 0, which therefore holds 
if and only if 

f' j 'f,i)h + rv* + v jk fjhh + 2^ k u Jt fu k = o , 

where we have replaced by Uji in the last term. But now 01]) shows 
that we can replace oJjif 1 with —fjiU 1 . This yields (1191) . Similarly, the 
equation 

= u l V l (f k f)f k + f k u J u k ) = --- 
shows that the final symmetry condition u l fi k cu k (uij—Uji) = reduces 
to (J2DD- □ 
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Corollary 4. Locally, a smooth function f with X < admits a smooth 
conjugate if and only if there is a smooth solution Ui of the equations 
©, ([T9D and PJJ. 

Proof. Symmetry of Uij is precisely the condition that Ui be exact and, 
therefore, locally of the form V^g for some smooth function g. □ 

Of course, we know that equations (jSD admit smooth solutions when 
X < so the only issue is whether we can find a solution for which 
(IT9|1 and (120]) are also satisfied. Also, if Ui is a solution then so is — w». 

4. Resolution of the equations: the generic case 

Throughout this section we shall suppose that X < 0. Recall that 
under this hypothesis / has four conjugate directions at each point, 
occurring in two pairs that differ only by sign. In other words, the 
solutions of the equations flSJ) have the form {±Wj, ±7^} for w f and rji 
smooth linearly independent 1-forms. Let us consider the expressions 

= f^ k hf J f k + P k hu J u k + 2Pf 3 k f l f k -2Pf J k u l u k 
p- = f* k fifth + f ijk fiVjVk + V ij fj k fifk - V ij fi k rim 

q - = p3kfJ. Vk + pk v . v . rjk+A p3f.kf. rjk 

According to Corollary H] and the discussion that immediately follows 
it, we now know that / admits a conjugate if and only if 

p + = q + = or P=Q= 0. 

These two possibilities are captured by the following theorem. 

Theorem 5. Locally, a smooth function f with X < admits a smooth 
conjugate if and only if 

p + p~ = o q + q~ = o (p + q) 2 + (pq + ) 2 = o. 

Proof. Evidently, the vanishing of these three quantities is equivalent 
to p ± = q ± =0. □ 

The condition p + p~ = was already resolved in [2]. We recapitulate 
and refine the argument as follows. Firstly, we write p + using normal 
coordinates (fl4|) to discover that 

(23) p + = p e + p bJ\LiJ2 

where 

(0A\ ^ e = f^f 333 + 2/3 2 (/l3 2 + f23 2 + ,/33 2 ) 

1 ' + (Miis - 2/ n 2 - 2/ 13 2 M 2 + (/s/223 - 2/ 22 2 - 2/ 23 2 V 2 2 
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and 

Po = 2/3/l23 - 4/13/23. 

In normal coordinates (771,772) = (±^i, T^)- It follows that 

(25) p - = Pe - p a u x u 2 
and hence that 

(26) = p 2 - p 2 u 2 u 2 2 . 

But, since Ui is subject to (j3J), we know that u;i 2 and a; 2 2 are determined 
in normal coordinates by f[T2"j) . In [2] we used this to eliminate Ui 2 and 
o; 2 2 from p e in (|24|) and then from p + p~ in (12"6"1) to discover by trial and 
error that Y 2 p + p~ could be written as an explicit Riemannian invariant 
in the derivatives of /, where Y is the invariant Z 2 — 2JX from (FTB]) . 
We can argue more systematically as follows. Firstly, we may obtain 
r]i from Ui without recourse to normal coordinates. 

Lemma 6. The conjugate direction rji is determined by the conjugate 
direction Ui via the formula 

(27) Wrn = 2f k f 3 uj k h + (Z- 2/ J ' fc /iA)wi - 2Jf i j u j . 

Proof. Since it is evidently coordinate-free, we may verify this formula 
in normal coordinates ( fl4l) . Substituting from (115j) we see that the 
right hand side of ( I2T1) becomes 

2(/l 3 / 3 Wl + /23/3W 2 )/ 4 + / 3 2 (/ll + / 22 )W 4 - 2f 3 2 f t j Uj 

In more detail, 



% 


right hand side of (1271) 


1 

2 
3 


/ 3 2 (/n + / 22 )wi - 2/ 3 2 (/ 11 u; 1 ) = / 3 2 (/ 22 - /n 
/ 3 2 (/n + / 22 )^ 2 - 2f 3 2 (f 22 u 2 ) = / 3 2 (/ n - / 22 
2(/i 3 / 3 ^i + / 23 / 3 c; 2 )/ 3 - 2f 3 2 (f 13 u 1 + f 23 u 2 ) -- 


)W2 

= 



On the other hand, from ffTBT) the left hand side of (|2"T|) becomes 

>//3 4 (/ll " / 22 ) 2 ^ 

and the whole of (1271) reduces to (7/1,7/2) = ±(wi, — w 2 ) depending on 
the sign chosen for the square root of Y. □ 

Note that since Y > when X < we could always insist of taking 
the positive square root of Y in (|27|) to obtain a consistent smooth 
choice of conjugate direction rji once is chosen. In any case, now let 
us consider p e in more detail. From (123]) and (|25|) we see that 

(28) Pe = i(p + +^)- 
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Note that p + does not see the sign of Ui and p~ does not see the 
sign of r]i. Moreover, interchanging Ui and 77, interchanges p + and 
p~. Hence, from ( 128|) we see that p e depends only on the derivatives 
of /. In principle, we could now use ( l2~Tj) to substitute for r\i in p~ . 
We conclude that Yp e is a polynomial in f i: fij, fij k , and Wj, which is 
actually independent of when fl5]) holds. Equation ([5]) may now be 
used to eliminate from Fp e leaving a polynomial in fi,fij,fijk- In 
practice, this is quite an intricate matter, which we consign to £TBl The 
result is: 

YPe = \Y{p + + p-) = \{ZS- 2XR + 2XY), 

where R and S are two further conformal invariants derived in £jA] 

Let us apply similar reasoning to some of the other quantities occur- 
ring above. From (I23p and ( 12 5 p we see that 

PoWiw 2 = \{p + -p~). 

As we have already observed, interchanging Ui and rji interchanges p + 
and p~ , hence changing the sign of p + — p~ . As is readily verified in 
normal coordinates, another quantity with this property is 

E = ^fiUjh'ut 

where e lJ is a choice of volume form, uniquely normalised up to sign by 
£^ k £ijk — 6. Specifically, if we further constrain our normal coordinates 
(fl4l) by requiring that e 123 = 1, then 

E = f 3 (J 22 - /n)wiw 2 . 

As above, it follows that we may use ([27)1 to eliminate 77, from 

YEpjj} x u 2 = \E{Yp + -Yp~). 

Moreover, this quantity is stable under interchange of Ui and 77^. It 
must be a polynomial in /j, fy, alone, which is given by: 

YEp u x u 2 = \E{Yp + - Yp~) = -\JXV, 

where this calculation is once more detailed in £jB] and V is one of our 
list of conformal invariants derived in £[A] But from Lemma [21 we have 
the identity 

(29) E 2 = -\J 2 X. 

We conclude that 

P = 8Y 2 p + p- = 2Y 2 (p + + p~ ) 2 - 2Y 2 (p + - p- ) 2 

= 2(ZS -2XR + 2XY) 2 + XV 2 . 

The vanishing of P is then our fourth conformally invariant condition 
(in addition to the first three (JSJ)), obtained in [2], for the existence of 
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a conjugate in the generic case X < 0. We now proceed similarly to 
obtain the two other conditions to provide a necessary and sufficient 
set of conditions. 

First we observe that Q = Y \/Y q + q~ is conformally invariant, where 
we use Lemma [6] to define rji by a choice of square root for Y. Cer- 
tainly it is a Riemannian invariant and we shall compute it in normal 
coordinates (I14p . According to the proof of Lemma 0, we may take 

v 7 ! 7 = h 2 (f22 - hi) Vi = wi r)2 = -w 2 , 
in which case 

q + = qxu)\ + q 2 U2 and q~ = qitui - q 2 C02, 

where 

qi = /s 2 /i33 + Anwi 2 + 3/i2 2 w 2 2 + 4/ 3 /i3 (A1 + /33) 

<?2 = /3 /233 + /222W2 +3/112^1 + 4/3/23 (/22 + /33J 

so that 

Q = Fv/Fg+g- = / 3 6 (/ 22 - /n) 3 (gi W - g 2 2 ^2 2 ) 

from which ui 2 and w 2 2 may be eliminated with (I12j) . The result is a 
polynomial expression in / and its derivatives. In terms of the various 
conformal invariants developed in §[A] it turns out that 

Q = \JZB -\JU -\ZS 2 

+ X(XZ 3 - JX 2 Z + 6W + \JM - ^ZXR + ^RS 

-fN+ \ZA — — IZK + g T + |C7 - g JD), 

as may be verified in normal form (|T4|) . 

The final condition (p + q~) 2 + (p~q + ) 2 = can similarly be expressed 
in terms of conformal invariants; although we do not attempt to write 
down the expression, we discuss how this can be done in £JE1 

5. Special cases 

5.1. Functions with a unique conjugate direction. Suppose now 
that / is a function that admits a unique conjugate direction up to 
sign. By f|T8l) . this occurs when X = and Y > 0. We first prove an 
analogue of Theorem [3j 

Theorem 7. Let u>j be a conjugate direction determined by with 
X = and Y > 0. Then the tensor field Uy is symmetric in its indices 
if and only if 

(30) e^ k f l u 3 {Jf k lm f l u m -2f kl f\f mn f m u n )) = 

(31) t %3k f l u 3 (jf k lm uj l uj m + f k l f l {r n f m f n + Z)) = 
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Proof. As in the proof of Theorem [3J is symmetric in its indices if 
and only if u % v^{uJij—Uji) = 0, where u l and are linearly independent 
vector fields. However, since X = 0, by Lemma [2j the vector field f^ojj 
is a linear combination of /' and u l and we have to use an alternative. 
A judicious choice turns out to be the vector field 

A short calculation using the identity 

shows that ey&i/ f^u k = —Z/J, which is non-zero by hypothesis (since 
Y = Z 2 ). In particular, v % has a non-zero component orthogonal to 
f l and to 1 . In order to bring this vector field into play, rather than 
differentiate the third equation from (J5J), we differentiate the equation: 

(32) e« fc /<^/fc'w { = 0. 
This gives 

(33) e ijk fiUjfkimU l + e ijk fmiUjfk l ui ~ m^v* = . 

First transvect this with f m . Then the resulting symmetry condition 
f m v % {oj m i — uJim) = holds if and only if 

But from (jlj), the last term can be replaced by uj m f im v % which is equal 
to [e lmn fif m r frUn/{f s fs)\f i:i fiUi (since e^ k f 3 u k f im u m = by iD). On 
multiplying through by J, we obtain the equation 
(34) 

Jj jk fi«>jfk lm fiUn ~ Je ijk fi l fifj m oo m u k - U" k J^ J h l J))J L,,^, = . 
However, from (132]) we deduce the identity 

Jf jm u m + (f M fkfi)^ - (f'hoJi)^ = . 

Indeed, the left-hand side is both orthogonal and colinear to the span 
of fj and LUj. On replacing J f jm u m by {f kl fk^i)fj - (f kl fkfi)ujj in the 
middle term of (1341) . we obtain (1301) . Similarly, on transvecting (1331) 
with u m , we conclude that the symmetry condition u m v l {uj m i— Ui m ) = 
is equivalent to (|311) . □ 

As for the generic case, we can summarise the conditions that / 
admits a conjugate as follows. 

Corollary 8. Locally, a smooth function f with X = admits a smooth 
conjugate if and only if there is a smooth solution of the equations 
©, d20D and flSTJ. 
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We can express these conditions in terms of the derivatives of / ei- 
ther by using invariant arguments, or by expressing them in normal 
coordinates. To do this invariantly, the following lemma can be em- 
ployed to eliminate quadratic terms in to 1 . 

Lemma 9. Suppose X = and Y ^ 0. Let Q 1 ^ be any symmetric 
form. Then 

(35) ZCpuVj = -ZCfUfi + ■IJO' 1 /,//".!).. + J\f k k Ql - Q kl f kl ) . 

Proof. Recall that E = e ljk fiUjfk e ooe satisfies E 2 = —J 2 X/2, so that 
(36) 

X = & E = & Jfjk^ + if'hfi^j-if'f^f^O, 

where the latter equality occurs since the LHS is both orthogonal and 
colinear to the span of fj and Uj. We then apply this to the identity 
given by transvecting f[iUJjfk k Qi] 1 = with f l ujK An alternative proof 
is simply to check that the formula holds in the Riemannian normali- 
sation. □ 

Equation (1301) can now be written in the form Q^UiUj = 0, where 

Qij = -e ikl f k (Jfr j f m -2flmf m (f nj fn)) 

-e ]kl f k (Jfrfm-2flmf m (rfn)), 

which, by Lemma [9] can be written as an invariant expression in the 
derivatives of /. However, it is more direct and somewhat simpler to 
just write out (130]) in the Riemannian normalisation. 

From the proof of Lemma [21 we see that X = implies that the 
product WiCo> 2 = 0. Thus (!30|) becomes: 

f 3 3 (co 1 2 -u 2 2 )(f 3 f 123 -2f 13 f 23 ) = 0. 

But since 7^0 (/ 2 2 — fn 0)> u i 2 ~~ UJ 2 2 = J and this is equivalent 
to 

/3 5 (/3/l23 — 2/13/23) = , 

which we recognize to be a multiple of V (which is given in normal 
coordinates by 4J 2 / 3 (/ 22 - /ii)(/ 3 /i 23 - 2/13/23))- Thus © and (J30D 
correspond to the conformally invariant condition V = 0. 

We give an invariant treatment of (13T|) as follows. Differentiate the 
right-hand identity of (1561) : 

= V,(J/ ife o; fe + (f kl fkfi)oJ j -(f kl f k ui)f j ) 

= 2(fuf l )fj k U k + J fijkUJ k + J fj k (jJik + fiklf k f l Uj + 2f kl fikflWj 

(37) + f M fkfi^ij — fikif k u l fj — f kl fik^ifj — f kl fk^ufj — f M fk^ifij ■ 
Note that for the moment we do not assume symmetry of Wy. 
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Recall the fundamental identities: u^Uj = f %3 fj and oj^fj = —f v u)j. 
Transvect (|37|) with w J to obtain: 

= -2f il ff jk f j fk + Jf tJ k^u k + J{ffJ)u ik 

• •/./•,/,•/./• , 7 , + 'I-! f h 'f,i.:(, + (f'hMhn ~ ./• / '• , ./•/.-•/./•,;-•' • 

From (E6D, J/jV = (f lm fi0o m )f k ~ (f lm fif m p k , so that 

J(/,-V> ife = {f lm f l U m )u ik f k -(f lm f l f m )u ik V k 

— —(f lm flUm)fikU k — (f lm flfm)fikf k , 

which gives the identity: 

From this, we deduce that (l3~Tj) has the equivalent form: 

e«' fc /i"i {-Jfk lm fif m - 2Jf M f lm f m + h l h{V mn fmfn + Z))=0^ 

(38) ^/.Wi {-(Tk + J( JV*(A/) - |A/V fc J)) = , 

where is the conformally invariant 1-form given by Theorem [T7] of 
Appendix lAl Even though JVfc(A/) — |A/Vfc J is not itself conformally 
invariant, its component orthogonal to the span of /j and Ui is, so 
the left-hand side of (1381) is conformally invariant. Now square this 
and use Lemma [9] to eliminate quadratic terms in Wj. We obtain an 
identity involving only the derivatives of /, which we identify in terms 
of conformal invariants as: 

(39) ffiV + §G + |F + iT - gZK - |ZA = . 

Corollary 10. Locally, a smooth function f with X = admits a 
smooth conjugate if and only if V = and ( 1391) are satisfied. 



5.2. Functions that admit infinitely many conjugates. When X 
and Y both vanish, the function / admits infinitely many conjugate 
directions. The following gives a complete description. 

Theorem 11. Suppose f is a smooth real-valued non-constant function 
such that its invariants X and Y both vanish. Then, up to scale and 
conformal transformation, f is one of the following 

(40) X\ log(xi + Xi + £3 ) arctan ' 



X2I Xi 2 + x 2 2 + x 3 2 



Proof. From ( 149]) we deduce immediately that = 0. But 



jy . . 1 1 1 / 1 (.'vrimi/'if f i'-W"'/"i i-Tii/ 1 t \>i rt , li 



the symmetric trace- free part of J 2 Vi\J 1 fj 
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whose vanishing is precisely saying that J~Vj is a conformal Killing 
field Vj all of which can be written down explicitly. Following [8], 

where Sj and j-j are arbitrary vectors, A is a arbitrary constant, and 
rrtij is an arbitrary skew matrix. We may invert 

Vj = U k fk)- l f 3 <=► /i = (v*^)" 1 ^ 

and inquire whether fj is closed. As a condition on 1/,-, this reads 
(41) VVfcVpVfl + 2V^V i] V fe = 0, 

the consequences of which are best viewed using a normal form for Vj 
such as those provided by Theorem [23 in §C\ Specifically, matrices of 
the form (IBTij) provide conformal Killing fields of the form 

. / d d d \ ( d d \ 

A PIq ^ X 2^ ^3^— + fJ,[X 2 ^ X 3^— 

V ox i 0x2 0x3/ V <yx 3 cte 2 / 

in accordance with the conventions of [8J. However, only when /1 = 
or A = is (1411) satisfied. When both vanish, we obtain the linear func- 
tions which are equivalent under scaling and conformal transformation 
to the first of (j4"0l) . Otherwise we obtain the second two, respectively. 
Matrices from the next group provide nothing new but matrices of the 
form ( ISTj) correspond to the conformal Killing fields 

/ d \ d d d 

li\x 2 - sc 3 -— - (xi 2 - x 2 2 - x 3 2 )- 2XiX 2 - 2XiX 3 -— 

V OX3 OX2 / OX\ OX2 ox 3 

and ( l4"Tj) is satisfied precisely when /i = 0. This gives rise to the final 
possibility for / in the list PUj) . □ 

In fact, all of the functions with X = Y = admit, not only infinitely 
many conjugate directions, but infinitely many conjugates. According 
to Theorem [TTJ it suffices to check this for the four cases (14"U|) . The 
first three of these are discussed in detail elsewhere in this article, 
specifically in §6.11 §6.3[ and §6.21 respectively. Finally, the functions 

x\ x 2 cos 6 + x 3 sin 6 

X\ 2 + X 2 2 + X 3 2 Xi 2 + x 2 2 + x 3 2 

form a conjugate pair for any 9. 
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5.3. Functions of two variables that admit a conjugate in M 3 . 

Let / = f(x2,Xs) be a function of two variables only. Then many 
conformal invariants simplify and in the case of a unique conjugate 
direction, the equations have a simple interpretation. As a first obser- 
vation, it is easily checked that X factors as a product: 

X = (A/)(/*ViJ), 

so that we also have 

Z = \fViJ + JAf , Y= (i/*ViJ- JAf) 2 . 

Furthermore, by its expression in the Riemannian normalisation, one 
sees that V = 0. In particular, the fourth condition for a conjugate: 
P = simplifies to 

ZS - 2XR + 2XY = . 

Now suppose X = and Y > 0. Then either Af = 0, in which case 
uj = (0, — fs, f 2 ) is, up to sign, the unique integrable conjugate direction 
and we are in the case of a planar function with planar conjugate, or 
/ l VjJ = and Af ^ 0. We can now exploit Theorem [71 Since (130]) 
is equivalent to V = 0, this is vacuous. However, ( )3T|) now comes into 
play. By going into the Riemannian normalisation, one sees that the 
third order terms of this equation vanish, and it becomes: 

(e^ k f i u J f k l f l )(f mn f m f n + Z) = 0. 

However, since Af 7^ 0, it is also the case that f mn f m fn + Z ^ and 
the equation becomes 

Let us write this out explicitly in coordinates: 

-wi/2/37t+wi/ 3 / 2 7i = o. 

But uj\ must be non-zero otherwise we are once more in the situation 
of a planar function with a planar conjugate whence Af = 0, contrary 
to our hypothesis. On combining this with the condition /*VjJ = 0, 
we obtain the simultaneous equations in / 2 ' fk and fe k fk'- 

{fsf2 k fk — f2fz k fk = 
fif2 k fk + hh fk = 

Since f 2 2 + fz 2 7^ 0, these only admit the solution f 2 k fk = h k fk = 0. 
But this implies that 

Vi(/*/fc) = ^ 1 1 V/| I = constant. 
The unique conjugate direction is thus given up to sign by 

^ = (v / /2 2 + / 3 2 ,0,0). 
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Furthermore this case occurs precisely when / satisfies the eikonal equa- 
tion ||V/|| 2 = constant. This should be compared with the example of 
a function having spherical symmetry as discussed in §6.3l below. where 
now the conjugate must satisfy an eikonal equation, even though there 
is no conformal transformation which sends concentric spheres to par- 
allel planes. 



In general, it is not the case that a function will admit a conjugate, 
even locally. For example, the function / = X1X2X3 has the property 
that X = 6f 2 . In particular X cannot be < on any open set, so that 
/ does not admit a conjugate on any open set. 

Recall from the Introduction that the pair (/, g) of a function and 
its conjugate define a semi-conformal mapping into M 2 . In the analytic 
category, such mappings arise (i) as the extension to the boundary at 
infinity of a harmonic morphism on the associated heaven space of the 
domain, see [5]; (ii) from local CR hypersurfaces in the standard Levi- 
indefinite hyperquadric in CP3, see [3]. The latter perspective leads to 
an explicit construction of semiconformal mappings from a holomorphic 
function of two complex variables, which, in a first form was given in 
[T2] then refined in [3J. In what follows, we highlight some particular 
cases of interest when a function / admits a conjugate function. 

6.1. Linear and quadratic functions. Any linear function / admits 
infinitely many conjugate functions, also linear; indeed the two invari- 
ants X and Y both vanish identically. The only quadratic function that 
admits a conjugate is, up to isometries and scaling, / = x± 2 — x 2 2 — x 3 2 . 
Note that / has an isolated critical point at the origin, however its 
conjugate g = xw/x^ + x 3 2 , although of class C 1 at the origin, is not 
smooth there. It is unknown if a pair of smooth conjugate functions 
(/, g) can have an isolated critical point. When they are harmonic and 
so determine a harmonic morphism, this is impossible [6] . 

6.2. Cylindrical symmetry. Let r 2 = x 2 2 + ^3 2 an d suppose that 
/ = / (r) so that its level sets are concentric cylinders. Then by solving 
the equations (j5), we obtain the conjugate direction:- 



whose four-valuedness corresponds to taking different signs for the 
square roots. Then for any branch, du — if and only if 



6. Some examples 




f 2 + vr 



c. 
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where C is a constant which is > 0. This has as first integral:- 
(42) f» = ^ + C, 



r 



where A > is a constant, and u is now given by 



^•2 ' i^«2 



A In | ^^" ^ } + VA + Cr 2 (C>0) 
VAlnr (C = 0) 



Then X = 2Cf'f"/r = -2AC/r 4 is < with the inequality strict 
provided neither of A nor C vanish. 

In fact we can integrate ()42p explicitly to obtain 

/ = 

The conjugate function is given by g = VCx\ — \J A arctan^/a^), 
interpolating between the two special case given by A = (/ = y/Cr) 
and C — (/ = v^41nr). In fact the mapping (/, p) has fibres which 
are helices lying on the cylinders r = constant. When C = these 
helices become circles lying in planes orthogonal to the xi-axis and 
when A = they become lines parallel to the ar-axis. Geometrically, 
we can interpret the four-valuedness of u as corresponding to the choice 
of a right-hand screw or a left-hand screw for the helices, together with 
a choice of orientation. In the special cases we obtain just two equal 
and opposite directions. 

6.3. Spherical symmetry. Let r 2 = X\ 2 + X2 2 +X3 2 and suppose that 
/ = f( r ) depends on the radial coordinate only. Then 

X = 2f(r) 2 (f"(r) + ££2 



so that if / is to admit a conjugate, the necessary condition X < 
forces / to be either constant or to satisfy the differential equation 

/" (r) + ™ = . 
r 

This has general solution / = Alogr + B, where A and B are arbitrary 
constants. For convenience, we take / = logr. Note that spherical 
symmetry implies that Y = and so there are infinitely many conjugate 
directions. In fact any conjugate function g must satisfy dg/dr = 
and ||Vg|| = 1/r. Thus g is determined by its values on say the sphere 
r = 1, where it must satisfy the equation \\Vg\ \ = 1. Such an equation 
is know as an eikonal equation and solutions are determined by initial 
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data on a hypersurface (i.e. a curve) in the sphere S 2 . It should be 
noted that the sphere S 2 does not admit a nowhere vanishing vector 
field and since we require 1 1 Vg 1 1 = 1 , then g cannot be globally defined 
on S 2 . Thus even though the function / defined on M 3 \ {0} admits 
infinitely many different conjugate functions in a neighbourhood of any 
point of its domain, the domain of any of these conjugate functions 
cannot coincide with that of /. 

6.4. An Ansatz. The following Ansatz provides a method of obtain- 
ing many pairs of conjugate functions. Let h(x, y) satisfy the partial 
differential equation: 

(43) (— X + 4y(— V + 4/i— = 

\dx J \dy J dy 

Then the functions 

f / = x 2 h(x 1 ,x 2 2 + x 3 2 ) 
\ 9 = x 3 h(x ll x 2 2 + x 3 2 ) 

are conjugate. For example, by taking h = (x 2 /y) + 1, we obtain 
the pair of conjugate functions of the Introduction. A straightforward 
calculation shows that the only product solutions h(x,y) = u(x)v(y) 
to (j4"3"j) . have the form 

be cx e y/l-cly 

h = . , 

1 + \[Y^y 

where b and c are constants. In fact, with reference to §5.1[ every 
solution obtained by this Ansatz satisfies X = 0. 

7. Invariants of the conjugate 

For a function / which admits a conjugate g, we can ask which of 
its properties are shared by its conjugate. More specifically, can we 
express the conformal invariants of g in terms of those of / ? For the 
invariant X, this turns out to be simply done. In order to be clear 
on which invariants are being considered, in this section we shall write 
X(f) and X(g) and so on, for the invariants of the respective functions. 

Theorem 12. If f admits a conjugate function g, then X(f) = X(g). 

Proof. In addition to flSJ), we have the identities: 

9' J J) ■ f^J, and g 'g, . 



X(g) = -2f kj g k fi gi + -{f kl f k9l ) 2 + -(/ w / fc /,) 2 
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Set Vi = €ijkf-'g k - Then we can decompose gry in terms of a symmetric 
basis: 

= j^(g kl fkfi)fifj + ^{g kl 9k9i)9i9j + -^{g^kVi^vj 

2 2 2 

+-p(g kl fkgi)f{ i g j ) + -j^{g fkVi)/^ + -^{g kl gkVi)g {i v j) 

1 2 2 

= -j(g h k)(J5ij - fifj - - -jfk(ig k fj) + -jfk(if h gj) 

1 2 

+-p{f kl fkgi){fif 3 - g l g 3 ) - -p{f kl fkfi)f(ig 3 ) ■ 

As a first application of this formula, we deduce the identity: 
(44) p gij -(f\)(g^) = 0. 

Furthermore 

9*9* = (g k k) 2 +^fk 3 g k f 3 gi+jfk 3 f k f h fi-j^(f kl fkgi) 2 -j^(f^ , 

which implies that 

In normal coordinates, on applying f fT2|) . the RHS equals 

-2g x 2 f 2 lx - 2g 2 2 f 22 2 + 2f 3 2 f 33 2 = 2/ 3 2 (/n + /33X/22 + /ss) , 

which is precisely X(f). □ 

Corollary 13. If f admits a conjugate function g, then for any eGR, 
the function f + eg admits g — ef as a conjugate and X(f + eg) = 
(l + e 2 ) 2 X(f). 

Proof. That / + eg and g — ef are conjugates, is easily checked. Then 

X(f + eg) = X(f) 

+eWhi + '.'/''//)./' 7 //■• - 2J(<?% - (fiWM 
+e 2 {4f i j f j g ik g k + \f, J g,g' k h + 2 gi j f j9 ik f k + 2fJg j f ik g k 

-APh + g^gij-ifif-Wj) 2 ]} 
■'''■{ HL'g, + \g, J f,)g"-gk - 2J(g ; 'f, J - (/',)(//',))} 
+e 4 X(g) . 

But the coefficients of the odd powers of e vanish on account of (jSJ) and 
( )44l) . so from Theorem [T2l we obtain 

X(f + eg) = X{f) + e 2 (X(f) + X{g)) + e 4 X(g) = (1 + e 2 ) 2 X{f) . 

□ 
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Note that if we view the pair (/, g) of a function and its conjugate as 
defining a semiconformal map into M 2 , then the replacement of (/, g) 
by (/ + eg,g — ef) amounts to multiplication of / + ig by 1 — it when 
we identify R 2 with the complex plane C. Indeed, semiconformality 
is preserved under conformal transformations of both the domain and 
co domain. 

To calculate the invariant Z(g) in terms of invariants of / turns out 
to be more challenging. In fact Z(g) depends on the choice of conju- 
gate direction, so that, in the generic case, the appropriate quantity to 
consider is the product \/Y Z (ui) Z (rj) . This can be calculated by the 
methods of §|B] to produce an expression involving third order deriva- 
tive of / which we don't attempt to write down. On the other hand, 
information about Z(g) can be obtained as in the above Corollary. 

Lemma 14. If f admits a conjugate g, then we have 

Z(f + eg) = (l + e 2 )(Z(f) + eZ(g)). 

Furthermore, 

Z(g) = ±Z(f + eg)\ e=0 ; 

that is, Z(g) = Zf(g) where Zf is the linearisation of the operator Z 
at f. 

In fact the latter part of the lemma is easily deduced directly from 

Z{g) = g lJ g igj + (g'giWj) = Phg 3 + (f l fi)(g j 3 ) , 

where, for a given / with V/ non-zero, the RHS is now a linear operator 
on g, which, since the principal term is the Laplacian, is elliptic. 

Proof. We have: 

Z{f + eg) = Z{f) + e(g i if i f j +2fig i f j + JAg) 

e 2 (2g^g i f j + f^g i9j + JAf) + e 3 Z(g) 
= Z(f) + e(Pf i9j + JAg) + <H.r'f,fj + JAf) + e 3 Z(g) 
= (l + e 2 )(Z(f) + eZ(g)). 

The last part of the lemma now follows from this formula, or as indi- 
cated above, directly from (jSJ). □ 

An interesting problem is to characterize those conjugate pairs that 
are 3-harmonic, i.e. conjugate pairs (/, g) satisfying Z(f) = Z(g) = 0, 
for then the mapping (/, g) determines a 3-harmonic morphism [TT]. If 
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both X and Y vanish, then so does Z and we have a complete descrip- 
tion in this case given by Theorem [TT1 Up to conformal transformation, 
the different conjugate 3-harmonic pairs are given by 

(xi,x 2 ), (|log(xi 2 + x 2 2 + x 3 2 ),arctan(x 3 /a;2)), 

, X\ X2 \ 

X\ 2 + X2 2 + X3 2 ' X\ 2 + X2 2 + X3 2 

More generally, by the homogeneity of Z(f) in /, the function / is 3- 
harmonic if and only if it satisfies the linearisation of Z at /: Zf(f) = 0, 
so that by Lemma H4"l Zf(f) = Zf{g) = is a necessary and sufficient 
condition for a conjugate pair (/, g) to be 3-harmonic. 

Appendix A. Conformal invariants 

Suppose / is a smooth function defined on an open subset U C M 3 . 
As usual, we denote the partial derivatives of / by subscripts 

fi = ?L A _ 92 f f k _ 93 f 

1 dx^ 13 dx i dx^ %3 dx i dxWx k ' 

Equivalently, we may regard these quantities as tensors obtained by 
repeated application of the flat connection Vj corresponding to the flat 
metric gr^. Suppose Q is a smooth non- vanishing function defined on 
U such that = fl 2 gij is also flat. If we let Tj = Vjlogfi, then it is 
well-known [6] that these functions are precisely the solutions of 

Vy y y 1„ yky 
i J- j J- i J- j 2 

and that all solutions are obtained by the conformal transformations 
of the round sphere S 3 viewed as flat-to-flat conformal rescalings via 
stereographic projection. Let Vj denote the metric connection for §ij 
and write 

/ = /, /i = Vi/, .[,, Y.Yjf. .... 
A conformal differential invariant of / of weight w is a polynomial 

I I (.9 3 i f \ fii fiji fijki ■ ■ •) 

in the derivatives of / and the inverse metric g lJ with the property that 
it is invariant under arbitrary coordinate transformation and 

(45) /(//''. /. J,, f, r f ijk , ...) = iVli^. f, /,, f, r f ijk , . . .) 

for all flat-to-flat conformal rescalings Q. As detailed in [9], this notion 
of invariance is the same as requiring equivariance under the action of 
SO(4, 1) on the 3-sphere, with M 3 <^-> 5" 3 by stereographic projection. 
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It is straightforward to write down explicit formulae for the effect of 
flat-to-flat rescalings on derivatives 

fi fi 

/ 4g N fij = fij ~ 2T (i/i) + 9ij^ k fk 

fijk = fijk — 6Y{ifjk) + 3g(ijY p f k ) p 

: + 6T(jTj/ fe ) - 3g(ijT k )T p f p - §T p Tpg (i .,/ fc ) 

with a view to verifying (|45p by direct calculation. It is difficult to 
find conformal invariants from this direct point of view. Certainly 
J = g^fifj = f l fi is an invariant of weight —2. Perhaps the simplest 
second order invariant is 

Z = ffifj + Jfj. 

It has weight —4 but it is usual to omit the powers of Q in verifying 
invariance (this is easily made precise by regarding the invariant as 
taking its values in an appropriate line-bundle). Specifically, as a linear 
combination of complete contractions it is manifestly invariant under 
coordinate transformation and 

f ij fjj = f'Lfr-T'.DpJ) /'•'/,/• JVf, 

Jf j j = JP, + Jr h f k 

whence 

f ij fjj + Jfj = rfifj + J?* 

as required. The familiar quantity 

(47) x = Wf 3 f k fk - ffif jk f jk + f'HFjf 

is a conformal invariant of weight —6. That it is a polynomial in the 
derivatives of / invariant under coordinate change is already manifest. 
Its conformal invariance, however, is most easily seen from the identity 
of Lemma r5J- 

JX + 12T ijk T^ k = 0, where T ijk = f^f^J. 

Here, recall that Uj is any solution of the equations §5§. We make take 
u)j = u>i to obtain a solution of the conformally transformed equations 
resulting from ( 146)) . Then 

fkiCj 1 = fkiuj 1 - rWfc + T'/^fc 

and so Tij k is a conformally invariant tensor (of weight —2). Notice, 
however, that Tij k is not an expression solely in / and its derivatives 
but also involves Uj. It may also be imaginary- valued. It is only in the 
combination Tij k T^ k that Uj can be eliminated using the relations (JSJ). 
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Of course, it is also possible to check the conformal invariance of X 
directly from the expression (j47j) . 

In the remainder of this section we construct an extensive menagerie 
of conformal differential invariants of /. It is possible, in principle [9], 
to list all such invariants. In practise, however, it is easier to construct 
invariants by a number of tricks (see [H]). Apart from the particular 
invariant V constructed below, these will turn out to be sufficient for 
our purposes. The new connection V« is related to Vj by 

V;0j = Vj0j - Tj0j - Tj<pi + gijT k (f) k 

when acting on an arbitrary 1-form <ftj. It follows that 

v*[jrVi] = n^vVi, 

which we will more conveniently express by saying if 0j has conformal 
weight —1, then 0; h-> V*0i is conformally invariant. Similarly, 

<t>j ^ V (i 0j) - lV k (p k g ij 

is conformally invariant when <pj has weight 2. Where J does not vanish 
we may consider the smooth 1-form J 1 / 2 /j. It has weight —1 whence 

is conformally invariant (of weight —4). As written here, this is not a 
polynomial but if we expand it we obtain 

\[V j J]fi + -JV'f, = pifih + ff.Pj , 

which is a perfectly good polynomial. It follows that this is an invariant 
whether or not J vanishes. It is our previous invariant Z. Another 
viewpoint on this construction is that f^Vj J + 2JV-,-/'- 7 is a conformally 
invariant bilinear differential pairing between /j and J. There are many 
such pairings on IR 3 as follows. 

Lemma 15. The following pairings are conformally invariant. 

S v ' 

weight v + w 
(V + lWV 4 - W0V;# 

11 v ' 

scalar of weight v + w — 2 
V 1p[iV + W<j}V[Hj}j] 

11 v ' 

skew of weight v + w 

(i;-2)[^V j) </)-|^ fc V fc 0] 
- iu0[V(iV>j) - ^ijV^fc] 

V v ' 

symmetric trace-free of weight v + w 



-0x0 i — y 

weight v weight w 
weight v weight w 

ditto i — y 

ditto i — y 
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Proof. These are all easily verified by direct calculation. Alternatively, 
we may employ evident variations on the trick used so far. For example, 
for non-vanishing ip and we may write the first pairing as 

(t)- v+1 i) w+l Vi[(t) v ^- w ], 

which is clearly invariant since cj) v il)~ w has weight zero. All of these 
pairings are similarly based on well-known conformally invariant linear 
differential operators. □ 

Notice that the bundles occurring in these pairings are irreducible in 
the sense that they are associated to irreducible representations of the 
orthogonal group. These are the bundles between which it is relatively 
straightforward to find invariant pairings. Here are two more examples 
that we shall need. 

Lemma 16. The following pairings are conformally invariant for if) of 
weight v and <pij being symmetric trace-free and of weight w. 

ip x (fiij i — y vipV^ij — (w + l)<j)ijS7 % i> 

" v ' 

weight v + w — 2 

-(w- 4)[0 fo -V fe) ^ - §0(ij0k)jVty] 

V v ' 

symmetric trace-free of weight v + w 

Proof. Easily verified by direct calculation. □ 

In fact, all the invariant pairings that we shall need may be con- 
structed from invariant linear differential operators. (There are, how- 
ever, many invariant pairings that do not arise in this way.) We are 
now able to list the almost all the conformal invariants that we shall 
use. 

Theorem 17. The following are conformal differential invariants of a 
smooth function f locally defined on M. 

J EE fU Z EE /-./•;./• + Jf{, 

x = ■if,\f J f ■'• ]),. - jp k f Jk + JU^ f 

If we now define 

Oi ee JViZ - 2ZViJ n ee JViX - 3XVi J 

4>ij = Jfij ~ ^f(ifj) k fk — \ jf k k9ij + ^f H fkfWij, 

then the following are also conformal invariants. 

R EE fGi S EE p Ti A EE ffV; B EE rVj 

D ee aWi T ee ftioiOj U = 
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// we now define 

Pijk = J^(i4>jk) ~ 30 (ii V fc ) J - lg(ijV l (t)k)i + lg(ij<Pk)N l J 

\j EE 2JV>y - V*J, 

then the following are also conformal invariants. 

F = p^^Xu G = (p ij XiXj K ee a'A, 

M EE r*Ai iV EE tTipV'Vifc W EE P ^ k Pij l (j) kl . 

Proof. We have already observed that J, Z, and X are conformally 
invariant. The remaining invariants in this theorem are manufactured 
from these basic ones by using Lemmata [TBI and [TBI as appropriate. □ 

There is one more invariant that we shall need and its construc- 
tion is slightly different. Let Q l i be any symmetric form and set 
v = e jkl (Jf k i Q lj - fQijfkmf^fi- Then the following identity holds: 

(48) Y&UiWj - Q lj ViVj) = ^Ev 

(recall that E = e ijk fiUjfk'ut). In the case when = p k f k -2f ik f k \ 
one may check that v is conformally invariant. It is convenient and 
consistent with [2] to define the related conformal invariant V = 4Jv. 
It has a different character to our previous invariants in that it changes 
sign under change of orientation. It is said to be an odd invariant. 

It is useful to record the conformal weight and homogeneity in / for 
each of the invariants of Theorem [T7| together with V:- 





J 


Z 


X 


R 


S 


V 


A 


B 


D 


weight 


-2 


-4 


— 


3 


-8 


-10 


-11 


-14 


-18 


-16 


degree 


2 


3 


4 




6 


7 


8 


10 


12 


11 




T 


U 




F 


G 


K 


M 


N 


W 


weight 


-18 


-22 




-18 


-18 


-14 


-16 


-18 


-18 


degree 


13 


15 




13 


13 


10 


11 


13 


13 



Any polynomial combination with consistent total weight will also be 
invariant. For example, the quantity Y = Z 2 — 2JX introduced in 
( Tl6l) is a conformal invariant of weight —8 (and homogeneity 6). Other 
evident invariants are not necessarily new. For example, it is easily 
verified by direct computation that 

(49) (jfifcj = \Z 2 - JX. 

This gives yet another verification that X is conformally invariant. 
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Appendix B. Invariant derivation of certain equations 

Our aim is to eliminate Ui from polynomial expressions of the form 
F(fi, fij, fijk, ....,u>i), given that the equations (jSJ) hold. We suppose 
that X < 0, so that in particular Y > 0. Recall that 

m = ^= {2(f kl fkVi)fi + (Jfk k - f kl fkfi)uJi - 2JfM} , 



gives the other conjugate direction, where an ambiguity of sign arises 
with the choice of square root. 

Lemma 18. Let be any symmetric form. Then 

(50) Y(Q ij UiUj + Q ij mm) = 2Q ij fif 3 (JX - Z 2 ) 

+2J 2 Q/{Zf l l - X) - -IJ-ZQ'-'f,, + UZQvtfhfj 

(51) VYQ^co^ = -Z&UU + ^Q"Lf/h + J\fu k Qi l - Q M fki) 

Proof. Both formulae can be deduced by skew-symmetrising over the 
indices of an appropriate 4-tensor. For example, to derive the second, 
consider the four tensor: T^i = fiUJjf k k Qi and apply the identity: 
Tujki] — 0. On transvecting first with then with up and applying 
(j3J), the result follows. □ 

Now let us find invariant proofs of some of the identities of §U Recall 

P + = p k hu 3 u k + fWfifjh - 2/^ w + ■>r<f/-f,f l , 

p- ee f^f iVjVk + p k fJ j f k -2f i Jf j %r ]k + 2fvf j k f i f k 

q+ EE fiKuiUjUk + fWfifjUt + lfVtffiUk 

q- ee f^rjirjjrjk + f ijh fifjVk + ^Pf k fWk ■ 
Theorem 19. The following identities hold: 

(52) Y(p + +p~) = ZS -2XR + 2XY 

(53) Y(p + -p~) = EV/J. 

where X, Y, Z, R, S, V are the standard conformal invariants and where 
E ee e ijk fi00jf k e uji. 

Proof. The first identity is an application of fl50l . where we have set 
Qij _ f l i k f k — 2f tk fy? . For the second, we apply (|48l) with symmetric 
form Qv = p k f k -2 f ik f k l. ' □ 

Note that both the LHS and the RHS of equation (1531) change sign 
under the interchange of the conjugate directions, the equation itself 
being well-defined and independent of this operation. 

The condition p + p~ = of Theorem now follows since 

AY 2 p + p- = Y 2 (p + +p-) 2 - Y 2 (p + - p-) 2 . 
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On applying ( |29|) . this gives the necessary condition P = of 
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2(ZS - 2XR + 2XY) 2 + XV 2 







Now consider the remaining conditions. We claim that we can use 
( 150]) and ( l5Tj) to write q l i k ujiUjUJ k as a linear form in Ui, where q^ k is 
any symmetric tensor. 

For this, first set Q l i = q l ^ k u k . Then from (15 Op . 

Y(t ih UiUjU} k = -Yq^WjCOk + 2q^ k f i f j co k (JX - Z 2 ) 



(54) 



+ 2J 2 q 3 i k uj k (Zf l l - X) - 2J 2 Zq» k f ij u k + UZ^tfftfw 



We now have to calculate Yq 13 r)ir)jU) k . For this we set Q lJ = x/Yq 13 rj k 
and apply (]5T]1: 

Yq %jk rjirjjU k = y/Y Q^ui^j 

= -Z^ h fifjV h + 2Jq ; ' k f l f/f ! r, + J\f/ gi lk v k - q ijk f ijVk ) 

where v % = y/Yrn = 2(f kl f k co l )f l + (Jf k k - f kl f k fi)oo t - 2Jf l k u k . On 
expanding the right-hand side and substituting into (154]) . we obtain: 

YqV k u i cj j u k = uAq^fjM-Y - 2Zf lm f l f m ) 



+ JqVj 
+J(Z 



Y + ZU mn fmfn)-2U mn fmfnf 



2f lm fif m )(2q^f/f l f k - Jq» k f jk ) - 2JZq> M fjf k f? 
2Jq jM f 3 fkfi m f m + J 2 Wq k M fi ~ J 2 q lkl f 3k fi - Zq^fjM 



+U 2 qi kl f j f k ™f m f l i + 2J\f/q l }\lY - qWfjJf) 

as claimed. 

We can now express Yq + by setting q^ k 



/ y and then adding 



Yq A 



Ui 



2Zf^f j f k (f lm f l f m ) 



Y + Z{r n fmfn) ~ 2{f mn f m fnf 

2f lm fif m )(2f ijk f/fifk ~ Jf jk f 3 k) ~ 2JZq^ l f j f k f, 



+Jf 1 

+J{Z- 

-2.r.r„ [>jp kl f J j),,frj,, l + J 2 f/f k kl fi - J 2 p kl f jk fi - zp kX i 3 hh 

+4J 2 P kl f j f k m f m f l i + 2J 3 (#/ fc w /, < - P kl f jk fi i ) + IP' JVJ) 



This has the form Yq + = a l uj^ where each a 1 is an explicit Riemannian 
invariant polynomial expression in fi,fij,fij k , which at each point is 
defined up to addition of an arbitrary linear combination: 



af + b[(f kl f k f l )u l + jp k u k }. 
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By symmetry, we must also have Yq = otr\i. Then the fifth condition 
YyY q + q~ = has the form r^uiirjj = 0, where r iJ is the symmetric 



form r u = (l/v7)aW. We can now apply (15T1) to give an invariant 
derivation of the quantity Q = Y 'VY q + q~ of §HJ 

The final equation of Theorem |5] is (p + q~~) 2 + {pq + ) 2 = 0. But 

H(p + q-) 2 + {P~q + ?} = {(P + + p-f + (P + ~ p-?}{{q + ? + (<T) 2 } 
-2(p++p-)(p+-^){(g+) 2 -(g-) 2 }, 

which we can see as a product of conformally invariant terms that we 
can deal with. First, multiply the whole expression by Y 3 \/Y . Then 
Y(p + + p~) is given by (1521) . whilst Y{p + — p~) is given by (155|) . On 
the other hand, 



YVY((q + Y + (q- 



which can be expressed using (l50l) above, whereas 



YVY((q+) 2 - (q 



^3 



can be expressed using (I48p . Note that the result involves E 2 , which 
by (129 p can be written in terms of conformal invariants of 



Appendix C. Normalising conformal Killing fields 

The conformal Killing fields on M. 3 form a finite-dimensional vector 
space on which 0(4, 1) acts via the conformal automorphisms of S 3 . It 
is the adjoint representation o(4, 1) and so the question of normalising 
a conformal Killing field up to conformal transformations comes down 
to finding canonical representatives for the orbits of this action. This 
is a question of linear algebra, which may be stated more generally 
as follows. Suppose we are given a real symmetric n x n matrix H 
of Lorentzian signature meaning that there is a real invertible n x n 
matrix such that 



(55) 



A l HA 



10 
0X0 
10 
-1 



Suppose iV is a real skew n x n matrix. We would like to find a real 
invertible n x n matrix A such that A t HA and A t NA are placed in 
some canonical form. For example, we may insist on ( |55|) for A t HA 
but following [8j |9] we normally prefer (written in block form) 



(56) 



A*HA 



1 

Id 

1 
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2!) 



where Id is the (n — 2) x (n — 2) identity matrix. 



Lemma 20. Suppose H is a real symmetric nxn matrix of Lorentzian 
signature and N is a real skew nxn matrix. Suppose that, regarded as a 
complex matrix, H~ l N has only one eigenvector up to scale. Then, the 
eigenvalue is zero, it must be that n = 3, and we can find an invertible 
real 3x3 matrix A such that 



A t HA 



1 

1 

1 



and A l NA 




2 
0-2 



Proof. Notice that 
H H> A 1 HA and 



N I— >■ A l NA 



H- l N I— >■ A~ l H~ l NA. 



Therefore, without loss of generality, we may suppose that H~ 1 N is 
in Jordan canonical form. Our hypothesis says that there is just one 
Jordan block with the eigenvalue A down the diagonal. But 

t^H^N) = tr(iy*(H*) _1 ) = -triNH- 1 ) = -t^H^N) 

so A = 0. In particular, the eigenspace is the same as the kernel of N. 
Suppose u is a non-zero vector in this kernel and consider 



u 



= {v s.t. u l Hv = 0}. 



Since 



rfHH^Nv = u l Nv = J^u = -v l Nu = 0, 



it follows that H~ l N preserves u 1 . The hypothesis that H^N has 
only one eigenvector up to scale now forces u G u L . In other words u is 
null, i.e. u l Hu = 0. It is well-known that 0(n — 1, 1) acts transitively 
on the null vectors. Therefore we may suppose that 



H 



1 
Id 




It follows that 



N 









and u 





M 



1 





where M is a skew (n — 2) x (n — 2) matrix. Therefore, 









M -r 




and (H^Nf 



r l M 
M 2 




-Mr 
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From the Jordan canonical form of H 1 N we see that, not only does 
its trace vanish, but also the traces of its higher powers. In particular, 

= ti((H~ l N) 2 ) = tr(M 2 ) 

and since M is skew it follows that M = and hence that rankiV = 2. 
Since the kernel of N is supposedly 1-dimensional, n = 3 is forced and 



Finally, if we take 



then 





' 





N = 





— r 




r 







" u- 1 


" 


A = 


1 










A* 



A l HA -- 
and so we can insist that fir 



1 
1 
10 



and ANA 








fir 





-fir 




—2 if we so wish. 



□ 



Lemma 21. Suppose H is a real symmetric nxn matrix of Lorentzian 
signature and N is a real skew nxn matrix. Then the eigenvalues of 
H~ 1 N lie on the real or imaginary axes. 

Proof. Suppose that x + iy is an eigenvalue, i.e. 

(57) H~ l N{u + iv) — (x + iy)(u + iv) for some u + iv ^ 0. 
Writing out the real and imaginary parts separately gives 

(58) H~ l Nu = xu — yv and H~ 1 Nv = yu + xv. 

We argue by contradiction, supposing that neither x nor y vanishes. In 



this case we see from ( 1581) that neither u nor v vanishes. Because iV is 
skew, we see from fl58|) that 

u l Nu 
v t Nv ■ 








u l HH x Nu = xu l Hu — yrfHv, 
v t HH~ 1 Nv = yv t Hu + xv t Hv. 

Therefore 

xu l Hu = yu l Hv = yv t Hu = —xv l Hv. 
Since we are supposing that i ^ 0, we conclude that u l Hu 
Again using (158]) . we now find that 

= u l Nv + v l Nu = u l H- l HNv + v t H~ 1 HNu = 2yu t Hu + 2xu l Hv, 
whence 

= xu l Hu — yu l Hv and = yu l Hu + xu l Hv. 



-v l Hv. 
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Therefore (x 2 + y 2 )u t Hu = and so u l Hu = 0. Bearing in mind our 
assumption that y ^ 0, we have found two real vectors u and v with 

u ^ 0, v ^ 0, u^Hu = 0, = 0, = 0. 

For H of Lorentzian signature this forces v — tu for some t 6 I. 
Substituting back into (|57|) and taking out a factor of (1 + it) gives 

H~ l Nu — (x + 

and hence that y = 0, our required contradiction. □ 

Lemma 22. Suppose H is a real symmetric nxn matrix of Lorentzian 
signature and N is a real skew nxn matrix. Suppose that H~ 1 N has 
a non-zero real eigenvalue A. Then —A is also an eigenvalue and we 
can find an invertible real nxn matrix A such that (in block form) 



AHA 



1 
Id 
10 



and ANA 








A 
M 




where M is a skew (n — 2) x (n — 2) matrix. 



Proof. Certainly, we may arrange that A t HA is of the required form 
and we shall suppose, without loss of generality, that H is already 
normalised like this. Write H~ l Nu = \u for Then 



= u l Nu = u'HH-'Nu = \u l Hu 
so u l Hu = 0. Therefore, by a suitable A we may arrange 



u 



and this forces H N 



Bearing in the mind that is skew, this implies 



N 



■ ■ A 
• 

-A 



and then H^N 





A 



-A 
• 
• • A 



It follows that —A is an eigenvalue, say H~ l Nv = —Xv for some v ^ 
and, reasoning as above, v t Hv = 0. Since u and v are not proportional, 
we may scale them so that u l Hv = 1. Finally, if we arrange that 



1 






then H~ X N = 
This immediately implies that has the desired form. 



-A 
• 
• A 



□ 
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With these Lemmata on hand we are now in a position to establish 
a general canonical form. As already mentioned, we shall prefer (156]) 
for A t HA. When n = 2 there is almost nothing more to do:- 



A HA 



1 

1 



and A l NA 



A 
-A 



simply because N is skew. It remains to observe that we can change 
the sign of A using 

" 1 
1 



.4 



but that A 2 is well-defined because the characteristic polynomial 

det(irtV-tid) = t 2 -X 2 
is invariant. The first interesting case is n = 3. 

Theorem 23. Suppose H is a real symmetric 3x3 matrix of Lorentzian 



signature and N is a real skew 3x3 matrix, 
invertible real 3x3 matrix A such that 



Then we can find an 



A l HA 



1 

10 

1 



and A l NA is 



(59) 









A 











-A 









or 



1 

71 






A 







' 





" 


-A 





-A 


or 








2 





A 










-2 






Furthermore, these three possible canonical forms are distinct apart 
from changing the sign of A in the first two cases and the coincidence 
of the first two cases when A = 0. 

Proof. If H~ X N has only one eigenvector up to scale, then Lemma 1201 
applies and we obtain the third case of ( l59l) . Else, Lemma I2T1 implies 
that either all eigenvalues are real or they are iX, —iX, for some A ^ 0. 

Firstly, let us suppose they are all real. They could still all be zero 
in which case the kernel N is at least 2-dimensional. But the rank of 
a skew matrix is always even so then N = 0. Otherwise, if A 7^ is a 
real eigenvalue, then Lemma [22] gives the first of ( |59l) . 

When iX is an eigenvalue, then 

H~ 1 N(u + iv) = iX(u + iv) 



implies that 



H^Nu = -Xv and H~ x Nv = Xu. 
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It follows that 
= u l Nu = u*HH- x Nu = -Xv/Hv 

= u l Nv + v l Nu = u l HH- l Nv + v^H^Nu = \{u l Hu - v*Hv) 
and so if A 7^ 0, we conclude that 

i/Hu = v t Hv and u l Hv = 0. 
In this case, by a suitable A we may arrange 



u 



1 

71 



and v 



from which the second of (j59j) is forced. Interchanging u and v changes 
the sign of A. Otherwise, the distinctions between these canonical forms 
is clear from the Jordan canonical form of H~ 1 N and its characteristic 
polynomial. □ 

It is easy to generalise these canonical forms tonxn matrices. The 
only one we shall need is the 5x5 case and we state it here. 

Theorem 24. Suppose H is a real symmetric 5x5 matrix of Lorentzian 
signature and N is a real skew 5x5 matrix. Then we can find an 
invertible real 5x5 matrix A such that 



A l HA 



1 

10 

10 

1 

1 



and A l NA is 

















A " 





















well-defined up to 


(60) 











/'■ 





(A,//) H> (-A,/i) or 










-A* 








or (-A,-/i), 




-A 

















or 



X/y/2 

-A/V2 -X/y/2 

ix 

-n 

X/y/2 



well-defined up to 
(A,/i) ^ (-A,/i) or (A, 
or (—A, —fi) or (/x, A) 
or (— /i, A) or (p, —A) 
or (-/!, -A), 
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or 



(61) 







-2 





2 
0^0 

-At 





well-defined up to ji ^ —fi. 



Furthermore, these canonical forms are distinct except for the evident 
coincidence of the first two cases when A = 0. 
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